Symmetry-protected topological phases of alkaline-earth cold fermionic atoms in one dimension 
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We investigate the existence of symmetry-protected topological phases in one-dimensional alkaline-earth cold 
fermionic atoms with general half-integer nuclear spin I at half filling. Using complementary techniques, we 
show that SU(2) topological phases are stabilized where the SU(2) symmetry stems from the existence of a 
metastable excited state in alkaline-earth atoms. On top of these phases, we find the emergence of topological 
phases with enlarged SU(2/ + 1) symmetry which depend only on the nuclear spins degrees of freedom. The 
main physical properties of the latter phases are further studied using a matrix-product state approach. We 
find that these phases are symmetry-protected topological phases, with respect to inversion symmetry, when 
/ — 1/2, 5/2, 9/2.., which is directly relevant to ytterbium and strontium cold fermions. 

PACS numbers: 71.1().Pm, 75.10.Pq 



Quantum phases of matter with exotic orderings have at- 
tracted much interest over the years. Prominent examples 
are topological phases which do not break any symmetry and 
cannot be characterized by local order parameters. Though 
all one-dimensional (ID) gapful phases have short-range en- 
tanglement, a topological phase can still be defined in ID by 
the presence of a symmetry which protects the phase.' With- 
out the symmetry, the phase can be smoothly connected to 
a trivial gapful phase without any quantum phase trans ition.EI 
In this respect, a topological phase with a symmetry is then 
called a symmetry-protected topological (SPT) phase.' ^ The 
Haldane phase^ of the spin-1 Heisenberg chain displays strik- 
ing properties which makes it a paradigmatic example of a 
SPT phase in ID. There is a finite gap above a unique ground 
state with periodic boundary conditions while spin- 1/2 edge 
states are liberated when the chain is open."* This phase ex- 
hibits also non-local string long-range ordering. These prop- 
erties of the Haldane phase have been shown to be protected 
by the presence of three global symmetries: the dihedral group 
of TT rotations along the x, y, z axes, time-reversal and inver- 
sion symmetries.^ 

A recent breakthrough leads to the complete classification 
of SPT of ID spin systems.** It relies on the determination 
of the projective representations of the underlying symmetry 
group which can be obtained using group cohomology. In par- 
ticular, for spin systems with a high SU(A^) symmetry, dis- 
tinct topological phases are expected from this classification.^ 
A natural question is the physical realization of these SU(A^) 
topological phases starting from realistic fermionic systems. 

In this respect, alkaline-earth like fermionic ultracold atoms 
seem to be very promising since they are the best candi- 
dates for the experimental realization of exotic high-symmetry 
many-body physics.""' In this context, the main interest in 
those atoms stems from the presence of a ground state ^Sq 
("g") and of a metastable excited state "^Pq ("e") between 
which transitions are forbidden. Moreover, both states have 
zero electronic angular momentum, so that the nuclear spin 
/ is decoupled from the electronic spin. The nuclear spin- 



dependent variation of the scattering lengths is expected to be 
smaller than ^ 10^^ for the ground state and ^ 10^'^ for 
the e state, and this results in fermionic systems with an ex- 
tended SU(A^ = 2/ + 1) symmetry.^ Such gases of strontium 
(*^^Sr) atoms (/ = 9/2) or ytterbium (^^^Yb, "^yb) atoms 
(/ = 1/2, 5/2 respectiv ely) h ave been cooled down to reach 
the quantum degeneracy.R^H] 

A simple model describing the low-energy properties of 
alkaline-earth ultracold fermionic atoms loaded into a ID op- 
tical lattice reads as follows;''* 

ida i i 

+ jY^[[T^f + {Tyf]+J^Y.^Tn\ (1) 

i i 

where ^ denotes the fermion creation operator at the z*'' 
site with orbital index I — g,e, and nuclear-spin states a = 
1, • • • , A^. In Eq. |Tj), ni = J^ia 4a, i'^ia, i denotes the occu- 
pation number on the site i, and 2Ti = J^ima 4a, i^imCma, % 
is the orbital pseudospin operator {a being the Pauli ma- 
trices). On top of the standard U(l)c charge symmetry 
{cia,i e*^c/Q,i), model features an SU(A^) symmetry 
(c/a,i — >■ ^pUapcip^i, U being an SU(A^) matrix), and a 
U(l)o symmetry in the orbital space spanned by the g, e states 
(Cg.ea. i ^ e^*^Cg.e „ i). The continuous symmetry of model 
^ is thus U(l)c X U(l)o X SU(Af). If one keeps only the g 
state, model ^ coiTesponds to the U(N) = U(l)c x SU(A^) 
Hubbard model. Its gapped Mott-insulating phases are known 
to be spatially nonuniform for all commensurate fillings.'^ In 
this respect, no topological phases can be formed in this one- 
band model. 

In this Letter, we show, by means of complementary meth- 
ods, that the situation is radically different if one includes 
the metastable excited state (i.e., the e state) of alkaline-earth 
atoms. The interplay between orbital and nuclear-spin degrees 
of freedom in this system gives rise to the emergence of differ- 
ent ID insulating topological phases at half filling (A^ atoms 
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per site). Two different classes of topological phases are stabi- 
lized when N is even, i.e., / is a half-odd integer. We find Hal- 
dane phases with integer orbital pseudospin-A^/2 and topolog- 
ical phases with enlarged SU(A^) symmetry that we fully char- 
acterize. In particular, we investigate the topological protec- 
tion of these phases with respect to the inversion symmetry by 
means of a matrix-product state (MPS) approach. It is shown 
that the topological phases of aUcaline-earth cold fermions are 
SPT phases when N/2 is odd, i.e., / ^ 1/2, 5/2, 9/2.., which 
is directly relevant to ytterbium and strontium atoms. 

Strong-coupling approach. A strong-coupling analysis 
along special lines of model ([T]l is useful to shed light on the 
possible existence of exotic Mott-insulating phases. In this re- 
spect, we consider the special case J = where the U(l)o 
orbital symmetry is enlarged to SU(2)o. The energy-spectrum 
for the single-site problem at half-filling (/z = UN) is then 
labelled by two integers p, q: 



E{p,q) 
D{p,q) 



^{p + q){p + q-2N)+'^{p-q){p- 

N\{N + l)\{p - q + If 
{N - p)\{N + 1 - q)\{p + ly.qV 
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where D is the degeneracy, n = p + q \s the number of 
fermions on one site, and T = {p — q)/2 is the spin of the 
orbital pseudospin operator T^. When J = 0, the model 
is U(2A^) symmetric. For U > 0, the lowest energy states 
have n = N fermions and states with different T are de- 
generate. These states transform in the self-conjugate anti- 
symmetric representations of SU(2A^). The resulting phase 
in the strong-coupling regime is known to be a spin-Peierls 
(SP) phase with bond ordering which breaks spontaneously 
the translation symmetry.'* In contrast, for U < 0, the low- 
est energy states have n = or 2N fermions and a charge- 
density wave (CDW) phase is stabilized.'** A second interest- 
ing line is [/ = and J < where the lowest energy states 
(Kji are + 1 degenerate (p — N and q = 0) and Ti is a 
N/2 pseudospin operator At second order of perturbation 
theory in | J| ^ t, we find a pseudospin- A^/2 antiferromag- 

netic SU(2) Heisenberg chain: Hctf = Jo J2i ' ^i+i ^i'^^ 
J„ = 8t'^/N{2N + 1)1 J|. As it is well known, the physics of 
the latter model strongly depends on the parity of A^.^ When 
A^ is odd, the phase is gapless with one gapless bosonic mode 
while the phase is fully gapped (Haldane gap) in the even A^ 
case. We thus find the emergence of an SU(2) topological 
(Haldane) phases when A^ is even, where the SU(2) symme- 
try stems from the orbital sector defined by the two g, e states 
of alkaline-earth atoms. We dub this phase, Haldane-orbital 
phase (HO), since this phase is different from the Haldane 
phase for (nuclear) spin degrees of freedom (see Fig. [T]where 
a cartoon represents this phase for A^ = 2). These Haldane 
phases are known to be topologically protected by a symmetry 
only when N/2 is odd, i.e., 7 = 1/2,5/2,9/2 for the simplest 
values.* 

A last interesting fine is J = 2NU/{N + 2) > with even 
A^. The lowest energy states correspond then top ~ q = N/2 
(N fermions which are orbital singlets) and transform into the 
SU(A^) self-conjugate representations described by a Young 




HO phase 



FIG. 1: (Color online) For N = 2: the blue (respectively yel- 
low) atoms are in the g (respectively e) state. The Haldane phases 
for nuclear spin and orbital (HO) degrees of freedom correspond 
to a (pseudo)spin-l Heisenberg chain. These phases present a di- 
lute ordeP that can be seen here as an alternation of sites with to- 
tal (pseudo)spin component ±1, diluted with an arbitrary number of 
sites with (pseudo)spin component 0. 



tableau with two columns and N/2 rows. At second order of 
perturbation theory in \U\ ^ t, we find an SU(A^) Heisen- 
berg spin chain in such representation. For N = 2, the spin- 
1 Haldane phase is formed among the nuclear spins and, in 
this respect, is intrinsically different from the HO phase with 
N = 2 (see Fig. 1). When N > 2, the physical properties of 
the ground state of the latter SU(A^) magnet are not known. A 
confinement of spinons is expected from the general classifi- 
cation of Ref 19 and a non-degenerate gapful phase has been 
predicted in the large A^ limit. 

Valence-bond-solid (VBS) approach. To get a good insight 
into the properties of the ground state of this SU(A^) Heisen- 
berg chain, we construct a series of model ground states, the 
VBS states,^' whose parent Hamiltonian is close to the origi- 
nal one in question. We start from a pair of the self-conjugate 
representations (characterized by a Young tableau with one 
column and A^/2 rows) on each site and create maximally- 
entangled pairs between adjacent sites (see Fig. |2]|. Last, 
we obtain the model VBS state by projecting the tensor- 
product states on each site onto the desired self-conjugate 
representation.^^ We expHcitly constructed the matrix-product 
representation 



Ai(mi)A2(m2)' 

{nii} 



■ Ai{m,i) ■ ■ ■ |mi,m2, 



■•) 

(3) 

of such a state for A^ = 4 (with the dimension of the A- 
matrices being 6 and {rrii} labeling the 20 physical states at 
each site) to obtain 'spin-spin' correlations exponentially de- 
caying with coiTelation length l/ln5 « 0.6213. The parent 
Hamiltonian which supports the above VBS state as the exact 
ground state is given by 



n 



VBS 



E 



qA qA 



108^^' +^6- 

(4) 

where Sf denote the SU(4) spin operators in the 20- 
dimensional representation. We observe that model (j4]) is 
not very far from the original pure Heisenberg Hamiltonian 
(with spin exchange Js) obtained by the t/J7-expansion. This 
strongly suggests that an SU(4) topological phase is stabilized 
in the strong-coupling regime with the emergent edge states 
belonging to the 6-dimensional representation of SU(4). From 
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the structure of the edge states, one can easily see that the VBS 
state found above belongs to one of the N topological classes 
protected by SU(iV) symmetry.'' 

It is interesting to see if it is robust even in the absence 
of SU(A^) symmetry or not. Recently, it has been demon- 
strated that the even-fold degenerate structure in the entan- 
glement spectrum is a fingerprint of the topological Haldane 
phase protected by inversion symmetry.^ One may think that 
the = 4 VBS state obtained above represents a topolog- 
ically robust state since its six finite entanglement eigenval- 
ues are all degenerate. However, this degeneracy is accidental 
and the corresponding topological phase is protected only by 
high symmetries. In fact, if the phase is protected by the link- 
inversion symmetry, the unitary matrix Uj satisfyin^Sl 

A^{m) ^e''^'UM{m)Ui (5) 

should be anti-symmetric.^ For the = 4 VBS state, Uj^ = 
Ui and the state becomes trivial only in the presence of such 
an elementary symmetry as the link-inversion. Since the di- 
mension of the self-conjugate representation is rapidly grow- 
ing (1764 for SU(8)), it is not practical to explicitly construct 
the MPS for > 6. However, it can be shown that the sym- 
metry of Ui is determined solely by that of the maximally- 
entangled (singlet) pair; if the maximally-entangled pair is 
anti-symmetric with respect to the interchange of the two con- 
stituent states, Ui is anti-symmetric and the corresponding 
MPS remains topological even without high symmetries. By 
investigating the form of the maximally-entangled pair, we 
conclude that for N = 2(2n + 1) (n > 0) the VBS state 
is in the stable topological phase protected by inversion sym- 
metry and not for — An. In particular, as in the HO case, 
we expect that the gapped SU(A^) phase realized along the line 
J = 2NU/{N + 2) is SPT when / = 1/2, 5/2, 9/2,etc. 

ffl ffl ffl 




maxima I ly-entang led 
pair 




FIG. 2: (Color online) SU(A'^) VBS states are constructed out of a 
pair of self-conjugate representations at each site. Dashed lines de- 
note maximally-entangled pairs, (a) SU(4) with 20-dimensional rep- 
resentation and (b) SU(6) with 175-dimensional representation. 

Low-energy approach. The low-energy effective field the- 
ory of the lattice model ([T]) is derived by expressing the 
standard continuum limit of the lattice fermionic operators 
ciai in terms of 2A^ left- and right-moving Li^, Ria Dirac 
fermionsPQa^ ~ i?/ae''=^^ + Ljc.e-''^^^, kp = 7r/(2ao) 
being the Fermi momentum and x = iao (ao being the lattice 
spacing). The continuum Hamiltonian then takes the form 
of 2N Dirac fermions coupled with marginal four-fermions 
interactions. The low-energy properties of the latter model 



are determined by means of a one-loop renormalization group 
(RG) analysis. This analysis has been done in the N — 2 case 
in Refs |25|2^ For instance, four different fully gapped Mott- 
insulating phases have been found when J = J^. On top of 
the two-fold degenerate phases with SP- and CDW orderings, 
the pseudospin-1 HO and the spin-1 Haldane phases, which 
have been identified in the strong-coupling analysis, persist in 
the weak-coupling regime as well. 

The RG analysis in the general N > 2 case is much more 
involved. In some regions of the phase diagram, we find that 
the one-loop RG flow is attracted along two isotropic rays 
with SO(4Af) symmetry which is the maximal continuous 
symmetry achievable for 2N Dirac fermions. These highly- 
symmetric rays signal the emergence of the CDW and SP 
phases in the general A^ case.'^In sharp contrast to A^ = 2, the 
one-loop RG flow for A^ > 2 features a region with no sym- 
metry restoration in the infrared (IR) limit. There is a separa- 
tion of the energy scales and one of the perturbations, which 
depends only on the SU(A^) spin degrees of freedom, reaches 
the strong-coupling regime faster than the others. A spin gap 
As is thus formed among the nuclear spin degrees of freedom. 
Below the energy scale of the spin gap E ^ As, the dominant 
pait of the effective interacting Hamiltonian only involves the 
remaining charge and orbital degrees of freedom and, when 
J = Jz, reads as follows: 

H^t ~ A {Tvgf + M cos ( ^SttXc/A^^c) , (6) 

where <i>c is a Bose field which accounts for the U(I)c charge 
degrees of freedom, and Kc is its Luttinger parameterl^S The 
low-energy properties of the charge degrees of freedom are 
captured by the sine-Gordon model at = SttKc/N and we 
expect that the charge sector is gapped away since Kc < N. 
The IR properties of model (j6]) thus depend only the orbital 
degrees of freedom and are described by the SU(2)jv confor- 
mal field theory (CFT) perturbed by its spin-1 operator (Trg)^ 
(g being the S U(2) matrix) with scaling dimension 4/ ( A^ + 2 ) . 
In this respect, the resulting low-energy field theory is ex- 
actly that of the spin-A^/2 SU(2) Heisenberg chain derived 
by Affleck and Haldane using the non-Abelian bosonization 
approach.l^i! We thus deduce the emergence of a Haldane-gap 
phase when A^ is even, i.e. for general half-integer nuclear 
spins. The resulting Haldane phase is identified with the HO 
phase which have been found already by the strong-coupling 
approach and is a collective singlet state formed among the 
orbital degrees of freedom. In stark contrast, within the weak- 
coupling approach, we could see no evidence for a similar 
SU(A^) topological phase of the nuclear spins for A^ > 2. In 
particular, along the J = 2NU/{N + 2) > line (A^ be- 
ing even), a SP phase is found in the RG approach instead 
of an SU(A^) topological phase expected from the strong- 
coupling argument. Therefore, a quantum phase transition 
necessarily occurs at intermediate couplings and it is tempt- 
ing to conjecture that the quantum critical point is described 
by an SU(A^)2 CFT. The situation may be understood as fol- 
lows. First, by symmetry, one can determine the low-energy 
field theory which is valid in the vicinity of this SU(A^)2 quan- 
tum critical point; TicS ^ 'Hsu(A')2 + (.9 — 5c)|Tr Gp, G being 
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the SU(iV)2 primary field which belongs to the fundamental 
representation of SU(iV). Following Refi28i the nature of the 
phases of the latter model can be inferred from a semiclas- 
sical approach; when g < gc, one has (TrG) 7^ 0, and the 
ground state is two-fold degenerate as a consequence of bro- 
ken translation symmetry (G — > — G). This may be identi- 
fied with the SP phase found in the weak-coupling limit. On 
the strong-coupling side g > gc, the semiclassical analysis 
now gives an SU(A^) matrix with the constraint TrG = 0, and 
the phase is translationally invariant. The resulting effective 
field theory is known to be the Grassmannian sigma model on 
U{N)/[U{N/2) X U{N/2)] manifold with a 6* = 27r, topo- 
logical theta term, which is massive.^^^ The latter is known 
to be the semiclassical field theory of the SU(A^) Heisen- 
berg spin chain in self-conjugate representations with two 
columns.^- Therefore, we conclude that the SU(A^) topolog- 
ical phase, identified within the VBS approach based on the 
strong-coupling Hamiltonian, emerges for g > gc- 

DMRG calculations. A density-matrix renormalization 
group^** (DMRG) is clearly called for to shed light on the exis- 
tence of this quantum phase transition for moderate couplings. 
We have rewritten model ([T]i as a A^-leg Hubbard ladder with 
additional rung interactions in order to get more efficient sim- 
ulations. Typically, we had to keep up to 1600 or 3200 states 
for convergence with N = 2 and = 4 respectively, and to 
get a discarded weight below 10^^, and we use open boundary 
conditions. 






J=4U/3 , 


^ U=1.5 




— U=1.8 




^ U=3.0 




" U=4.0 
i i i t 





20 



30 



1.6 

S 1.4 
c 

■a 1-2 
-a 

g 0.8 

o 

0.6 



- (b) ' 


^ a=l 


1 




• ■ a=2 






^ a=3 






" a=4 






U=6 J=8 I \. 


1 


1 









10 



20 



30 



FIG. 3: (Color online) (a) Bond kinetic energy for various U > 
along the line J = 4(7/3 (N = 4) on L = 36 chain, the hopping t 
is fixed to one. (b) local densities rii^a (see text for definition) in the 
topological phase (U — 6 and J — 8) showing the existence of edge 
states. The dashed line is an exponential fit with a correlation length 
of around 3 lattice spacings. 



a uniform non-degenerate one. This is in agreement with our 
conjecture based on the low-energy analysis and the strong 
coupling one. In order to confirm that the non-degenerate 
phase corresponds to an SU(4) topological phase, we plot in 



Fig. gb) 



the local densities n, 



; ^/a, i'-'la, i 



for each 



flavor a — 1, . . . , 4 in this phase. Clearly, edge states are 
present with an exponential profile and a correlation length 
of the order of 3 lattice spacings for t/ = 6 and J = 8 (see 
fit). Moreover, DMRG data randomly shows different kinds of 
edge states, in the sense that for a given edge, two flavors out 
of four give the same local densities. The same phenomenon 
occurs in the Haldane phase of spin-1 chain where a simu- 
lation done in the 5*°* — sector will generate arbitrarily 
an edge state with Sz — ±0.5 at one boundary. Now, since 
there are 6 ways to choose two flavors out of four, our nu- 
merical data confirm that edge states are 6-fold degenerate 
in agreement with the VBS prediction that they belong to 6- 
dimensional SU(4) representation. We also observe that these 
edge states are absent in the SP phase. The precise location of 
the phase transition and its universality class require extensive 
large-scale simulations and are beyond the scope of the Letter 
Now, we turn to the illustration of the HO phase for N = 2 
since the case = 4 would be quite demanding numeri- 
cally.'^ In Fig. Qa), we show the presence of edge states in 
a finite chain of length L = 64, similarly to what is found in 
the Haldane phase of the spin-1 chain. In Fig. |4|b), we also 
plot pseudo-spin correlations (taken from the middle of the 
chain) that are short-ranged and almost identical to spin cor- 
relations measured in a S* = 1 Heisenberg chain. Overall, we 
confirm the existence of the HO phase in this region. 
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FIG. 4: (Color online) A = 2 DMRG numerical data for a 
chain of length L = 64 with [/ = and J/t = -12 (a) lo- 
cal average {Tz{x)) vs position x in the ground-state with total 
Tz — 1 showing evidence of localized edge states; (b) pseudospin 
correlations r^^(x) = (r^(i/2)r^(L/2 + x)) and Tprr,{x) = 
{T+{L/2)T-{L/2 + x))/2 in the ground-state with = exhibit 
an SU{2)o symmetry, are short-ranged and almost identical to spin 
correlations {Sz{L/2)Sz{L/2 + x)) measured in 5 = 1 Heisen- 
berg chain. Bulk and edge correlation lengths are close to 5 lattice 
spacings. 



In Fig. |3ja), we plot the bond kinetic energy for various 
U > using A^ = 4 along the special line J = 2NU/{N + 
2) — 4C//3 for a chain of length L = 36. As U increases, 
we clearly see a quantum phase transition from a SP phase to 
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